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Abstract: Dark matter trapped in the Sun produces a flux of all flavors of neutrinos, 
which then reach the Earth after propagating out of the Sun and oscillating from the 
production point to the detector. The typical signal which is looked at refers to the muon 
O • neutrino component and consists of a flux of up-going muons in a neutrino detector, the 

q^j major source of background being atmospheric neutrinos. We propose instead a novel 

signature, namely the possibility of looking at the tau neutrino component of the dark 
matter signal, which is almost background-free in the downward-going direction, since 
the tau neutrino amount in atmospheric neutrinos is negligible and in the down-going 
baseline atmospheric muon-neutrinos have no time to sizably oscillate. We analyze the 
prospects of studying the downward-going tau neutrinos from dark matter annihilation 
(or decay) in the Sun in Cherenkov detectors, by looking at hadronic showers produced in 
the charged-current tau neutrino interactions and subsequent tau decay. We discuss the 
various sources of background (namely the small tau neutrino component in atmospheric 
neutrinos, both from direct production and from oscillations; tau neutrinos from solar 
corona interactions; the galactic tau neutrino component) as well as sources of background 
due to misidentification of electron and muon events. We find that the downward-going tau 
neutrinos signal has potentially very good prospects for Mton scale Cherenkov detectors, 
the main limitation being the level of misidentification of non-tau events, which need to be 
kept at level of percent. Several tens of events per year (depending on the dark matter mass 
and annihilation/decay channel) are potentially collectible with a Mton scale detector, and 
a 5 a significance discovery is potentially reachable for dark matter masses in the range 
from 20 to 300 GeV with a few years of exposure on a Mton detector. 
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1 Introduction 

It is well established by means of astrophysical and cosmological probes that almost 23% 
of the matter present in the Universe is in the form of a non-luminous component, the 
so-called Dark Matter (DM). Although the evidence for DM in cosmic structures is of 
gravitational origin, it is expected that DM particles may produce a large variety of di- 
rect or indirect signals, which represent the subject of a large number of experimental and 
theoretical studies. While direct searches exploit the direct scattering of DM particles on 
the nuclei of a low-background detector, indirect signals look for the products of DM self- 
annihilation (or decay) occurring in the galactic halo, in the extragalactic environment or 
in those celestial bodies, like the Sun and the Earth, where they can be gravitationally 
trapped and accumulated before starting the annihilation (or decay) process. All annihila- 
tion (or decay) products are absorbed by the bodies, except neutrinos: a flux of high-energy 
neutrinos, in general composed by all three flavours, can then emerge from the Sun or from 
the core of the Earth as a signal of the presence of the trapped DM. 

The typical way of looking for a signal from DM annihilation in the Sun or in the Earth 
is to search for an excess of upward-going muons produced by charged current interactions 
of the i/« flux produced by DM annihilation. The advantage to look at the upward-going 
neutrino flux though its muon conversion is represented by the large conversion area offered 
by the rock below the detector, and by the well-established and efficient experimental 
techniques to detect high energy muons. The main source of background is therefore 
represented by the f^-component of atmospheric neutrinos. 
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A great deal of work has been devoted to these studies, see for instance Refs. [1], and 
Refs. [2] for applications to the case of neutralino DM. Limits on the flux of through- 
going muons from DM annihilation inside the Sun, the Earth and at the Galactic Center 
have been set using data from water Cherenkov detectors like Super-Kamiokande (SK) 
[3] and from neutrino telescopes like AMANDA and IceCube [4]. Competitive limits can 
also be obtained considering the stopping muons, which are usually the dominant signal 
for low-mass DM particles. This analysis has been carried out in the context of a light— 
neutralino signal at SK in Ref. [5] and in a more general framework in Ref. [6], where 
fully contained events were also considered. Recently, it has been discussed in Ref. [7] the 
ability of studying the low mass DM range using track events in the DeepCore array of the 
IceCube detector, that can reach an energy threshold as low as 10 GeV. Detectors designed 
for other physics searches, such as liquid argon and magnetized iron calorimeter detectors, 
are able to achieve a much better energy and directional resolution than water Cherenkov 
detectors and thus could be used for detecting neutrinos from DM annihilations inside the 
Sun [8]. Moreover, the possibility of detecting DM through electron neutrinos in liquid 
scintillation experiments like KamLAND has also been considered, see e.g. Ref. [9]. 

In this paper we instead intend to propose and study the feasibility of a different 
signature: downward-going tau neutrinos produced by DM annihilation in the Sun. The 
basic motivation is the following: while DM annihilation typically produces similar amounts 
of all neutrino flavors, atmospheric neutrinos are largely dominated by the v e and 
components, the u T 's being very suppressed. Neutrino oscillations transform a fraction of 
the original atmospheric Va into u T , but this phenomenon is basically inoperative in the 
down-going direction, due to the much larger oscillation length (of the order of the Earth 
diameter) as compared to the average production height of atmospheric neutrinos, that 
for definiteness we fixed to be 15 km in the vertical direction (see e.g. Ref. [10] for a 
more precise discussion of the production height). This fact implies that in the downward 
direction the v T produced inside the Sun by DM annihilation (or decay) represents an 
almost background-free signal. 

We wish to recall that the v T component in the DM neutrino signal coming from the 
Sun is basically unavoidable: i) the mechanism of production of DM annihilation is hardly 
flavor-sensitive to the level of producing only v e and v„\ ii) neutrino oscillations on the 
baseline of the Sun-Earth distance in any case average out any production of into a 
fraction of v T able to reach the Earth. This implies that the amount of the DM-signal v r J s 
reaching the Earth from the Sun are of the same order of magnitude as the DM-signal f^'s, 
but they are basically background-free. Therefore they may be a competitive discovery 
channel as compared to the standard up-going za/s signal, which is affected by a much 
larger background represented by atmospheric f«'s. 

The background sources for the v T channel are represented by the very small intrinsic 
v T production from cosmic-rays interactions in the atmosphere (atmospheric zVs) or by 
residual atmospheric Vn oscillation, negligible for down-going but more relevant for almost 
horizontal fluxes, by oscillation of v e and produced by cosmic-rays interactions in the 
solar corona (which represent an irreducible background for our DM signal, since it comes 
from the same source, the Sun) and by a galactic plane v T flux, mainly given by oscillation 
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of Vp_ produced in cosmic-ray interactions with the interstellar medium (this source of 
background is, in principle, reducible by angular cuts on the galactic plane position). We 
will show that these sources of backgrounds are, by themselves, not posing significant limits 
to the DM signal, which therefore represents a potentially viable novel possibility. 

We will instead show that the limiting factor for an analysis of the downward-going 
z/ T 's is the ability of a neutrino detector to identify the tau neutrinos. We will concentrate 
on water Cherenkov detectors, in order to be able to access the relatively low energies 
which are relevant for DM studies [E v < TeV) (in the case of DM with mass from few 
tens of GeV to few TeV) . Since next generation apparata of Mton scale are currently under 
study, our predictions are particularly suited for future water Cherenkov detectors like the 
Hyper-Kamiokande (HK) project [11]. 

A signal from tau neutrinos, specifically in connection with DM searches, have been 
also studied in Ref. [12], where the possibility to detect v T from gravitino DM decay in 
the halo of our Galaxy was discussed. High energetic neutrinos (E v > TeV) coming from 
astrophysical sources have been considered in Refs. [13, 14], where the flux of v r arising from 
neutrino oscillation was discussed. The ability to detect atmospheric tau neutrino events 
through an ultra-large liquid argon detector is described in Ref. [15]. An experimental 
measurement to detect tau atmospheric neutrinos by measuring the energy spectra of 
neutrino induced showers has been suggested in Ref. [16]. 

The scheme of the present paper is the following. In Sect. 2 we will discuss the 
neutrino fluxes coming from DM annihilation and the background fluxes: the atmospheric, 
solar corona and galactic neutrinos. The class of signals relevant to water Cherenkov 
detectors and a discussion about the experimental limits in detecting v T in those type of 
detectors is given in Sect. 3. The calculation of the hadronic contained events is then 
presented in Sect. 4, both for the signal and for the various sources of background. The 
capabilities of Mton-scale water Cherenkov detectors on the discovery of the v T signal from 
DM annihilation coming from the Sun is derived in Sect. 5, where we discuss also the level 
of sensitivity that can be achieved from this detection channel on the DM-proton scattering 
cross section a v . Conclusions are given in Sect. 6. 

2 Neutrino fluxes 

Let us start by discussing the relevant neutrino fluxes produced by DM annihilation in the 
Sun. The decay case is easily obtained along the same line of reasoning. We are interested 
not only on the source spectra, but also in the propagated fluxes which reach the Earth, 
and that go through both energy redistribution inside the solar medium and oscillation 
processes. We will discuss and treat all these phenomena as in Ref. [17]. 

In addition to the DM signal, we need to discuss the relevant sources of background 
for the v T signal, which arise from atmospheric, solar corona and galactic cosmic-rays 
interactions. 



-3- 



2.1 Dark Matter 

The neutrino flux at the detector, coming from DM annihilation inside the Sun and for 
each of the three neutrino flavours, is defined as: 

r m dNf 



dE v 

f 

where BRf is the DM branching ratio into the annihilation final-state channel /, dNl /dE u 
are the neutrino spectra for each channel /, T Q denotes the annihilation rate inside the 
Sun and d is the distance between the Earth and the Sun. For definiteness, in our analysis 
we will consider two specific annihilation channels: i) direct annihilation into neutrinos 
with flavor-blind branching ratios, i.e., BR Ue p e = BR^p^ = BR Ur p T = 1/3; ii) annihilation 
into rf leptons with BR Tf = 1. More general situations are easily implemented. The 
DM fluxes dNl/dE u are calculated as described in Refs. [17, 18], considering all physical 
processes that can occur to the neutrinos from the production point inside the Sun to 
the detector at the Earth: vacuum and matter neutrino oscillations, neutral current and 
charged current interactions on the Sun's medium. We have fixed the neutrino mixing 
angles as [19]: sin 2 #23 = 0.5, sin 2 612 = 0.304 and #13 = 0. The annihilation rate r is 
calculated as in Ref . [20] . For the DM velocity distribution in the rest frame of the Galaxy 
we have assumed a Maxwellian distribution with v = 220 km s _1 and we have fixed the 
local DM density po to be 0.3 GeV cm -3 . The annihilation rate depends on the DM-proton 
scattering cross section, which we will fix to benchmark values in our analysis, as specified 
in the following. 

Since we are interested in the (downward going) neutrino flux coming from the Sun, it 
can be important to know the amount of time spent by the Sun at different zenith angles 
above the horizon. This is relevant to determine the duty factor of the signal and could 
also be exploited for optimizing the signal-to-background ratio, in the case the detection 
of the v T flux could be correlated to the Sun's position in the sky. If we consider a detector 
located at a latitude ip, then the motion of the Sun is described by the following expression 
[21]: 

cos 9z = sine/? sin 5 + cosy? cos 5 cos#ha (2-2) 

where 9z is the zenith angle, 5 is the solar declination and #ha is the hour angle. The 
value of 5 can be obtained using an approximate expression, 

t = 23.45° sin (360°(^-80)N 

V 365 J V ' 

where N^ &y denotes the number of the day, starting on the first of January. The hour angle 
#HA is zero at the local solar noon and indicates the time that has passed since the Sun 
was at the local meridian. In the left panel of Fig. 1 we report the apparent motion of the 
Sun, considering a detector located at a latitude cp = 36°, that is roughly the latitude of 
the Kamioka site. In the right panel of Fig. 1 the amount of time (in hours) spent in one 
year by the Sun in each zenith-angle bin is shown, for the same latitude: this determines 
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Figure 1. Left panel: Apparent motion of the Sun as a function of time (hour and day of the 
year), for a detector located at the latitude ip = 36°. The regions denoted by different shades 
denote the value of the cosine of the zenith angle (shown by the numbers in labels) occupied by the 
Sun at a given day and hour. Positive values of cosOz refers to the Sun above the horizon. Right 
panel: Amount of time (in hours) spent by the Sun in each zenith-angle bin in one year, for a 
latitude (p = 36° (solid line). The dashed gray line shows the time spent by the galactic center in 
each zenith-angle bin. The dotted orange line represents the averaged duty factor of the galactic 
plane convoluted with the position of the Sun, see text for more details. 

the yearly duty-factor for the signal. Since z/ T 's produced in the galactic plane are one of 
our sources of background, in the same right panel of Fig. 1 we also show the time that 
the galactic plane spends in each zenith angle bin. We have averaged over 120 points in 
the galactic plane (b = 0) and for each bin above (below) the horizon we have imposed 
the condition that also the Sun position is above (below) the horizon, even if not in the 
same angular bin of the galactic plane. This is done to determine the duty-factor of the 
galactic-plane contribution correlated with the presence of the Sun in the upper (or lower) 
hemisphere. For this reason, the integral of the galactic plane duty-factor over the total 
zenith angle bins does not turn out to be one year. Just for definiteness, we also show in 
Fig. 1 the duty-factor of the galactic center without considering the position of the Sun. 
In this case the duty-factor is zero for cos 9z > 0.5. 

2.2 Backgrounds 

There are various sources of background for the v T signal from DM annihilation. First of 
all we have the v T atmospheric background coming from oscillation of atmospheric v^. This 
form of background is the dominant one for upward-going neutrinos, but it is extremely 
reduced for the downward-going case that we are going to consider. For the atmospheric 
neutrino, we have used the Honda fluxes [22], considering the new release from February 
2011, where the implementation of the JAM interaction model [23] for the low energy 
interactions of cosmic rays and air nuclei was implemented. We considered the azimuth- 
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angle-averaged flux at Kamioka site with mountain over the detector. 

A second form of background is given by the intrinsic v T contribution to the atmo- 
spheric flux coming from decay of charmed particles. This has been computed in Ref. [24] 
and can be described by the following parameterization: 



log 



10 



n3 d4u , ( GeV^ 



dti v \ cur s sr 



A + Bx-Cx 2 - Dx A , (2.4) 



which is valid for 10 2 GeV < E v < 10 6 GeV, with x = log 10 {E„ [GeV] ) , A = 6.69, B = 1.05, 
C = 0.150 and D = -0.00820. We have extrapolated the flux down to Ef™ = 3.5 GeV, 
that is the minimal neutrino energy required for the production of a tau lepton. In order 
to check our extrapolation, we have compared the fluxes that we obtain in this way with 
the ones reported in Ref. [25]. We found that for energies E u < 10 GeV our results are in 
well agreement with [25], since the v T fluxes are dominated by the oscillated atmospheric 
neutrinos in that range of energies. For 10 GeV < E v < 10 2 GeV, our results are equivalent 
to [25] for the horizontal direction, while they differ by at most a factor of two for the down- 
going direction (cos9 z = 1). We have checked, however, that the events coming from this 
energy range account for roughly the 11% of the total v T events, that we will present in 
Sect. 3. In the same Section, we will see that the actual experimental background is highly 
dominated by the misidentified events. Thus, we expect that our extrapolation of the 
intrinsic fluxes between 10 and 10 2 GeV does not have a relevant impact on our analysis 
(we will comment more precisely on this in Sect. 5). We want, moreover, to stress that the 
calculation of the intrinsic v T component for energies E v > 10 2 GeV suffers from sizable 
uncertainties. In fact, the atmospheric showering parameters are not precisely known [26] 
and in the literature different charm production models are present: quark gluon string 
model [27], recombination quark proton model [28] and perturbative QCD [29]. The fluxes 
that we use in this article are based on a perturbative QCD approach [24], but depending 
on the specific models, the predictions for the fluxes might change of up to one orders of 
magnitude. We refer to Refs. [26, 30, 31] for a detailed discussion on this topic and we 
will comment in Sect. 5 on possible implication of this uncertainties on DM searches. In 
the following, when we refer to 'atmospheric' v T flux we will always imply the sum of the 
oscillated fluxes from atmospheric Un and the intrinsic contribution. 

The third form of background is represented by the neutrino flux produced in the 
solar corona by cosmic-ray collisions. This has been studied in Ref. [32], where a suit- 
able parameterization of the flux of electron and muon neutrinos (j = v e ,v^) has been 
determined: 

^ ^o^^TU^Vcm 2 ^-^, (2.5) 



dE v u I + A{E y [Ge\}) 

which is valid for 10 2 GeV < E v < 10 6 GeV. For the numerical values of the coefficients 
Nq, A and 7 we refer to Ref. [32]. Also in this case we have extrapolated the neutrino 
fluxes down to £7™ m = 3.5 GeV. During their propagation to the Earth, the electron and 
muon neutrinos produced in the solar corona oscillate and generate u T . 
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Finally, the fourth possible source of background is given by the fluxes of tau neutrinos 
from the Galactic plane. A complete discussion of this topic is presented in Ref. [33, 34]. 
The v T flux can be parameterized as: 



rgi = 9 x l(T 6 (GeV cm 2 s sr^^GeV])" 2 - 64 , (2.6) 

which is valid in the energy range 1 GeV < E v < 10 3 GeV. As discussed in connection with 
Fig. 1, for a detector located at a latitude (p = 36° this contribution starts to be sizable for 
cos9z ^$ 0.5. In this range of zenith angles the atmospheric and intrinsic fluxes are by far 
the dominant ones and thus, for the v T signal coming from the Sun, the galactic neutrinos 
do not represent a troublesome source of background. We have nevertheless included this 
contribution in the calculation. 

Fig. 2 shows the total neutrino fluxes (integrated from E™ in = 3.5 GeV up to E™ n = 
10 4 GeV) as a function of the zenith angle cos 6z- It is well visible how the atmospheric 
v T background is sizebly reduced for cos 9z > as compared to the cos 9z < case 
and how the v T flux is much smaller than the v e and fluxes. This behavior is the 
main motivation for our proposal on downward-going tau events. In fact, as already 
discussed above, the DM v T signal from the Sun is comparable, in size, with the signal. 
Nevertheless, the up-going muon channel (measured through up-going muons) has a much 
larger atmospheric background than the down-going tau channel, as is clearly seen in Fig. 2. 
An up-going muon signal which occurs to be completely dominated by the atmospheric 
neutrino background, could instead dominate the atmospheric background in the down- 
going tau channel. In Fig. 2 we also show the solar corona and the galactic v r fluxes. We 
notice that solar corona neutrinos become non-negligible mainly for cos #z~0.5, while the 
galactic neutrinos are more important for angles with cos#z < 0.5. For the down-going v r 
signal, atmospheric neutrinos represent the main source of background, while the galactic 
and solar corona neutrinos give a sizable contribution for small and large zenith angles 
respectively. In both cases, they are very suppressed, therefore offering potential chances 
to a signal to emerge. 

A few examples of downward-going v T fluxes as a function of the neutrino energy 
E u are shown in Fig. 3, together with the total downward-going v T background fluxes: 
sum of atmospheric, both intrinsic and from oscillation, solar corona and galactic fluxes. 
To allow for a realistic comparison between the fluxes, in summing the galactic and solar 
contributions to the atmospheric ones, we have multiplied by the fraction of the time that 
the Sun or galactic plane spend in each zenith angle bin with respect to the total time 
that the Sun spends above the horizon in one year. We have integrated over zenith angles 
with co's.Qz > 0. This represents the angular range that we will consider throughout the 
paper. For the signal, we show the cases for DM mass of 10, 100, 1000 GeV, and for 
the DM spin-independent scattering cross section a p we have used a benchmark value of 
10~ 41 cm 2 (with a p we will denote throughout the paper the spin-independent scattering 
cross section, while with <jp D the spin dependent one). The left and right panels refer to the 
two benchmark annihilation channels: direct annihilation into neutrinos and annihilation 
into tau leptons. In the case of direct annihilation into neutrinos, we clearly notice the 
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Figure 2. Atmospheric neutrino fluxes as a function of the zenith angle cos6>z at Kamioka site. 
The dotted (blue), solid (red) and dashed (green) line refer to the v e , and v T fluxes at the 
detector (with oscillation included) . The solid horizontal (black) line denotes the v T flux from solar 
corona interaction, while the dotted horizontal (orange) line refers to the level of the galactic in- 
flux. Note that these two fluxes (differently from the atmospheric contributions) arc present only 
when the Sun or the galactic plane occupy each specific angular bin. Fluxes have been integrated 
from E™ in = 3.5 GeV. 



line a E v = m x together with the degraded-energy tail due to neutrino propagation in the 
Sun's interior. Also from Fig. 3 we can see that the signal is significantly dominant over 
the background in the down-going v T channel, especially at high energies where the tau 
cross section is not suppressed by mass threshold. In the same figure, we also show the 
flux of down-going atmospheric which (together with the down-going z/ e ) will produce 
a source of background events for the down-going u T signal in water Cherenkov detectors, 
due to misidentification of the muon and electron events in the detector, as will be discussed 
in Sect. 3. 



3 Signals at water Cherenkov detectors 

The experimental signals of high-energetic tau neutrinos have been analyzed in Ref. [35], 
where the double bang and lollipop signatures were considered. These signatures, however, 
are distinctive only for energies above the PeV range. For lower energies, between TeV 
and PeV range, it might be possible to tag taus which decay to muons, if the neutrino 
interaction vertex occurs within the detector [35]. 

For the range of energies we are interested in (GeV-TeV), charged-current v T interac- 
tions in water Cherenkov detectors will lead to multiple Cherenkov rings and the possibility 
of individually identifying these events is currently based only on statistical methods. In 
Ref. [36, 37] the SK Collaboration has employed neural network and maximum likeli- 
hood techniques to successfully discriminate tau neutrino events. They concentrate on the 
hadronic decays of tau leptons {BR^i ~ 64%), since they have a more spherical topology 
than the backgrounds. We will exploit the same type of signature in our calculations. The 
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Figure 3. Fluxes of v T arising from DM annihilation inside the Sun: the left panel refers to the case 
of direct annihilation into neutrinos with even branching ratios in the three flavours, while the right 
panel stands for annihilation into tau leptons. In the left panel, the continuum part of the spectrum 
is in units of GeV -1 m~ 2 s _1 , while the line is in units of m~ 2 s _1 . The DM spin-independent 
scattering cross section a v is fixed at the value 10 -41 cm 2 and the cases for DM masses of 10, 100 
and 1000 GeV are reported. The solid (red) line shows the downward-going flux of atmospheric 
muon neutrinos v^, while the long-dashed (green) line refers to the downward-going background 
flux of v T (sum of atmospheric, both intrinsic and from oscillation, solar corona and galactic fluxes). 
The fluxes have been integrated over zenith angles with cos8z > 0. 



primary backgrounds to v T charged-current events are neutral-current and charged-current 
events from u e and atmospheric neutrinos. A number of event selection criteria [36, 37] 
can be applied to reduce these backgrounds. Using a likelihood analysis and a neural 
network, the SK Collaboration achieved an efficiency (with respect to the total number of 
events in the fiducial volume) of 43.1% and 39.0% for tau events identification, respectively. 
The electron and muon events are misidentified as tau events with a percentage of 3.8% 
and 3.4% for the same two statistical analyses. Since the atmospheric v e and Vn fluxes 
are usually larger than or of the same order of the DM fluxes, as can be seen in Fig. 3, 
the misidentification has a relevant impact on the actual performance of water Cherenkov 
detectors in identifying v T events. 

In our analysis we have focused on the favourable situation in which the track events 
from Un charged current interactions will always be detected and correctly identified. For 
this reason, in the misidentified background events we will only considered neutral-current 
events from v e and atmospheric neutrinos and charged-current events from u e interac- 
tions. Note that in principle also v e and from solar corona interactions will contribute 
to this form of background, but they are negligible with respect to the atmospheric v e 
and Vp,. Another possible form of background events is represented by neutral-current 
events from v T (atmospheric, intrinsic and solar corona contributions) and by v T charged- 
current interaction, with tau decaying into electron. We will neglect these events since 
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Figure 4. Iso-contours of number of downward-going v T hadronic events in the plane of detector 
exposure (in kton x year) and DM mass m x . The left panel refers to DM annihilation into neutrinos, 
while the right panels stands for annihilation into tau leptons. The elastic scattering cross section 
on protons (relevant for capture in the Sun) is fixed at the benchmark value of 10 -41 cm 2 . The 
solid (red), dotted (blue) and dashed (black) lines refer to 5, 10 and 50 events in the detector, 
respectively. The two horizontal (pink) lines denote the exposures that can be reached by a 0.5 
Mton detector, like HK, in 1 (solid) and 10 (dotted) years. 

we will numerically see that the background is typically dominated by neutral current is e 
and i>u and by charged-current u e misidentified events. Finally, we wish to comment that 
also Vp, and v e fluxes from DM annihilations would contribute to the signal events, again 
through misidentification. We will not include this type of contribution, which is marginal 
in determining the amount of signal events. 

4 Contained hadronic events 

The class of signal events we are considering is represented by the charged-current pro- 
duction of a tau lepton, followed by its hadronic decay. The hadronic showers of the decay 
produce Cherenkov rings in the water detector. In this Section we discuss our determi- 
nation of the number of hadronic events, and the calculation of the class of events which 
contribute to the background through misidentification, namely background events from 
v e and Un neutral current interactions and from v e charged current interactions. 

4.1 Signal events from hadronic tau decay 

For the calculation of contained tau events, we will use the derivation of Ref. [38], where 
the concept of visible energy is introduced. On this topic, see also Refs. [14, 16, 39]. For 
this category of events, the visible energy E v { s is the sum of the energy of the broken 
nucleon and the hadronic energy Ef l 2 of the tau decay. The total number of contained 
events for charged current v T interactions is: 
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r^r r dree 

N^ c L B = M det N y x dE vis / dn 77(g) — /_ , (4.1) 

where M det is the detector mass, N y the number of years of exposure and 77(g) is the on- 
source duty factor. For the signal, we use the duty-factor discussed in Sect. 2.1 and shown 
in the right panel of Fig. 1. For the background events we have conservatively considered 
the case in which the events are classified only as upward/downward. Thus, the background 
is not filtered through the duty-factor of Fig. 1, instead it is always present when the Sun 
is above the horizon (77 = 1/2). In a more optimistic scenario in which the detector would 
be able to correlate the direction of the observed Cherenkov rings of the hadronic showers 
with the position of the Sun, a reduction of the background would be possible. We will 
nevertheless neglect this optimistic possibility here and we will consider only the more 
conservative case in which the events are identified as upward-going or downward-going. 
The function dl^ c /d£l dE vls in Eq .(4.1) is defined as: 



dT? 



dn dE v 



S,B 



dE v / dE T 



dttdE v 



T,^ C (E T , E v 



SB 



dE vl 



(4.2) 



where S and B denote signal and background, respectively, d(p UT / dQdE v is the v T flux 
coming from DM annihilation or the background flux of atmospheric, intrinsic and solar 
corona tau neutrinos. The function S^ c , defined as: 



Yf?(E T ,E v )=N A 



(4.3) 



quantifies the number of interactions, with N A being the Avogadro's number, Z and M 
the fraction of proton and neutrons present in the detector and daf,'™ jdE T the v T charged- 
current cross section on proton and neutron, for which we adopt Ref. [40], where the 
correction due to the finite tau mass is implemented. For the parton distribution functions 
(PDF) we have chosen the MSTW 2008 NNLO [41] varying Q 2 from 5 GeV 2 till 10 4 GeV 2 . 
For lower Q 2 we have frozen the PDF to the values assumed at Q 2 =5 GeV 2 . For a more 
refined calculation, it should be taken into account also the non-perturbative [42] evolution 
of the PDF in the low Q 2 region, Q 2 < 2 GeV 2 . Moreover, also target mass correction [43] 
could have some impact for very precise estimates. Nevertheless, we have neglected these 
two latter corrections, since they are beyond the precision required for our study. We have 
checked that with our assumption we achieve an agreement better than 10% with the total 
cross section reported in Ref. [40]. 

The function dTh/dE v - m in Eq. (4.2) is defined as the decay rate of tau leptons into 
hadrons, provided that the energy of the hadronic decays £^2 is equal to the tau lepton 
energy minus the final neutrino energy E Uf , with E v = E v — E v [ s : 



dE h 2 JT -, h 5(E hj2 



dE. 



I1.2 



(E T - (E u - Eyi,,))) G(E U - E v 



(4.4) 
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Figure 5. Number of downward-going v T hadronic events (solid lines) as a function of the DM 
mass m x , for a scattering cross section on proton a = 10™ 41 cm 2 and for a detector exposure 
Ma et N y =1 Mtonxyear. We show also as horizontal lines the neutral current events expected from 
the atmospheric u e (lower dotted blue line) and (solid red curve), the charged-current events 
from v e (upper dotted line) and the charged-current events expected from the background v T . 



The decay rate dnh/dE^^ can be obtained from the v T spectra produced in the decay: 

dn h 1 dn h E h2 . . 

with Zh = , (4.5) 



r 



dE h>2 E T dz h E. 
where 

dn h dn v E v - E vis 

- — = - — with z v = . (4.6) 

dzh dz v E T 

The expressions of dn v jdz v are given in Refs. [14, 24, 44]. Note that the function dn^jdE^i 
is normalized to the total branching ratio of tau into hadrons, BR^ = 0.64. 
Using the previous definitions, we can rewrite Eq. (4.2) as: 



dx c S 



dn dE v 



S,B 



dE u / dE T 



fT . J I?min 



X C T e (E T ,E v )±^, (4.7) 



with an analogous part for antineutrinos and with E™ m = max[m T , E v — E^ x ^\. We have 
considered a lower limit on the visible energy equal to the minimal neutrino energy for tau 
lepton production: = 3.5 GeV. A specific detector should choose the minimal visible 

energy able to maximize the number of correctly identified tau events with respect to the 
misidentified ones. This quantity could turn out to be different from our choice of E 1 ^ 1 , 
but this would be highly detector dependent and should be obtained through dedicated 
Monte Carlo analyses. We fix the upper limit on the visible energy to E™ x — 10 4 GeV for 
the atmospheric events, or to E™ x = m x for the DM events. 

The number of signal events expected for a detector of a given exposure (expressed 
in kton x year) is shown in Fig. 4 as a function of the DM mass, for the two cases of 
DM annihilation into neutrinos (left panel) and annihilation into tau leptons (right panel) . 
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The elastic scattering cross section on protons, relevant for capture in the Sun and that 
determines the size of the annihilation rate Tq, is fixed at the benchmark value of 10 -41 cm 2 . 
In all our analyses we are considering that equilibrium between capture and annihilation has 
been reached, as typically occurs for the Sun. The solid (red), dotted (blue) and dashed 
(black) lines refer to 5, 10 and 50 events in the water Cherenkov detector, respectively. 
The two horizontal (pink) lines denote the exposures that can be reached by a 0.5 Mton 
detector, like HK, in 1 (solid) and 10 (dotted) years. We notice that for a Mton-scale 
detector, the expected number of down-going v T signal events can reach the level of 50 or 
more, depending on the DM mass. For a prolonged exposure of a decade or more, more 
than a hundred of signal events is potentially under reach. For the annihilation channels 
under study, this signal is most sensitive to DM masses in the range from 30 GeV up to 200- 
300 GeV. For lighter DM, the signal is reduced, since the available energy to be transferred 
to the hadronic showers in the tau channel is reduced (we remind that only neutrinos with 
energy E v < m x are produced by the non-relativistic annihilation process). For heavier 
DM the signal fades away since capture is less efficient and moreover absorption processes 
of higher energy neutrinos in the Sun start to become operative [17, 18]. These properties 
can also be observes in Fig. 5, where the expected number of downward-going v T hadronic 
events are shown (as solid lines) as a function of the DM mass m x , for the same scattering 
cross section on proton of Fig. 4, for a detector exposure M& ct N y =1 Mton x year. In 
the same figure, we also show the charged-current events expected from the background v r 
(horizontal dashed line) . This figure demonstrates that the specific pv-induced background 
is negligible. 

4.2 Background events from v e and v„ neutral current interactions 

As we have already mentioned before, a non-negligible source of background is represented 
by the atmospheric v e and v„ neutral current interactions, since a fraction of the Cherenkov 
ring they produce are misidentified as hadronic tau events. The rate of these events is 
defined as: 



N^L = M det N y x / ^ dE vis dn r}{6) 
J Ej?i n J 

where dl^/dUdE^ [38] is: 



dU dE v 



(4.8) 



dCl dE v 



v dndE v 



^{E Vls ,E v ) + {u 



v). 



(4.9) 



The function is the analogous of Eq. (4.3) for neutral current interactions, for which 
we use the deep-inelastic cross section of Ref. [45]. Note that, for neutral current events, 
the fraction of energy that is transferred to the cascade is given by y = E Y - 1S /E V , with E v 
being the initial neutrino energy. 

The number of neutral current v e and v„ events for one Mton x year exposure are re- 
ported in Fig. 5, as horizontal lines. The lower dotted blue line shows the neutral current 
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events expected from the atmospheric v e , the solid red curve refers to z/„. Together with 
the charged-current v e events (discussed in the next Section) , these class of events are not 
directly comparable to the signal events, since they represent a backgound only when they 
are misidentified as tau hadronic events. Fig. 4 shows that these classes of events pose a 
problem if they are not controlled at a level better than a few percent. This level of misiden- 
tification is foreseeable, since SK reconstruction and analysis algorithms are already able to 
reduce the misidentification at the level of less than 10% [36, 37]. A misidentification level 
around 1%, together with an efficiency of reconstruction of hadronic tau-events increased of 
about 70-80% would bring the down-going v r signal to be a competitive search. A possible 
way to reduce the misidentification background could be obtained by considering specific 
tau hadronic decays. Indeed, the tau lepton decays mainly producing multiple pion events: 
about 40% of the hadronic decays are given by r~ — > 7r°7r~z^ T . This channel produces two 
electromagnetic cascades from ir° — > 77 decay and one muon track from tt~ — > decay. 
The main experimental challenge would be to find suitable cuts to statistically distinguish 
these events from the CC/NC u e and multiple pion productions. 

4.3 Background events from u e charged current interactions 

Another source of background for the hadronic tau decay is represented by the charged- 
current v e channel. We calculate background events due to z/ e charge current interactions 
as [38]: 



ATfl = M det N y x / dE vis I cM n{0) 



B .-aet^- , —v* , — dndR 



£f TOT (£ vis ) + ( u ^u), (4.10) 



where Eg ^ot is the analogous of Eq. (4.3) for the total charged-current cross section of 
u e , for which we used the deep-inelastic expressions of Ref. [45]. Note that, in this case, 
the visible energy E V [ S is equal to the full initial neutrino energy E v . 

The number of charged current v e events for one Mtonxyear exposure is reported again 
in Fig. 5. We notice that this is the category of events which (through misidentification) 
represents the largest source of background. 

5 Detectability and statistical significance 

To quantify the discovery reach of present and future water Cherenkov detectors, we use 
the statistical significance defined as the signal-to-noise ratio: 

<*wn- <"> 

See also Ref. [46] for a complete discussion on the statistical significance. We have studied 
the behaviour of q in two cases: an ideal case in which no misidentification is present and 
the detector efficiency for tau leptons is almost 100%, and a more realistic case in which 
both the misidentification and the detection efficiency for taus are considered. In the ideal 
case we have: 

Adcal = N? C \ B . (5.2) 
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In the presence of misidentification of the electron e™ ls and muon e™ ls events and detection 
efficiency e T for the detection of tau hadronic events, we instead have: 

Realistic = e T N c r % + ef* + ef s {N? C \ B + N?\ B ) . (5.3) 

Our results are reported in Fig. 6, 7 and 8. Fig. 6 shows the iso-contours of statis- 
tical significance for the detection of downward-going v T hadronic events in the plane of 
detector exposure (in kton x year) and DM mass m x . The left panel refers to DM an- 
nihilation into neutrinos, while the right panels stands for annihilation into tau leptons. 
The elastic scattering cross section on protons (relevant for capture in the Sun) is fixed at 
the benchmark value of 10 -41 cm 2 . The solid (red), dotted (blue) and dashed (black) lines 
refer to significance of 2, 3 and 5 a. The two horizontal (pink) lines denote the exposures 
that can be reached by a 0.5 Mton detector, like HK, in 1 (solid) and 10 (dotted) years. 
This plot refers to an ideal detector with full efficiency for the detection of tau hadronic 
events (e T = 100%) respect to the neutral and charged current electron and muon events 
( e ™ ls = e ™ s = 0)- We notice that in this ideal case a 5 a significance of signal detection is 
present for a wide DM mass range (from about 20 to 300 GeV) with 1 year of exposure of 
1 Mton detector. Indication at the 2-3 a level are possible up to 500-700 GeV for the case 
of direct annihilation into neutrinos, while it extends to larger masses (and drops al small 
DM masses) for the case of annihilation into tau leptons. Clearly, statistical significance 
scales with the scattering cross section a p . Note that a more precise analysis would require 
to take into account the theoretical uncertainties on the intrinsic v T fluxes [26, 30, 31]. The 
intrinsic contribution constitutes roughly 21% of the total v T events calculated as described 
in Sect. 4.1. For an intrinsic flux one order of magnitude higher than the one considered 
here [24], the curves for the statistical significance ? in Fig. 6 would move up by a factor of 
1.3 for the rf annihilation channel and for a DM mass m x between 10 and 20 GeV. The 
curves will remain roughly unchanged for higher masses and for the vu annihilation chan- 
nel. We will moreover see that in the realistic cases of our interest (with misidentification) 
the dependence on the intrinsic v T flux is highly reduced. 

A more realistic situation requires to take under proper consideration the efficiency for 
reconstruction of tau events and the misidentification of the v e and Un events. This is shown 
in Fig. 7, where we reduce the efficiency to e T = 40% and we allow a misidentification of 
4%, both in the electron and muon channels, not far from those already achieved by the 
SK analysis [36, 37]. We notice that now large statistical significance is recovered with a 
much larger exposure: nevertheless, a 10 yr period of data taking reproduces the same level 
of reach as the one discussed for the optimal case in Fig. 6. Fig. 8 shows the capabilities 
in the case of better performance: e T = 70% and e™ ls = e™ ls = 1%. In this case, a few 
years of exposure would suffice to cover almost the whole DM mass range (with noticeable 
differences for different annihilation channels) for our benchmark value of the scattering 
cross section a p . Note that also for e™ ls = e™ ls = 1% the background is still dominated 
by misidentified events, see Fig. 5, and the uncertainties on the v T intrinsic flux (which 
account for roughly 21% of the total number of events) does not alter our conclusions. 

The dependence of the statistical significance with the detection efficiencies is reported, 
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Figure 6. Iso-contours of statistical significance for the detection of downward-going v T hadronic 
events in the plane of detector exposure (in kton x year) and DM mass m x . The left panel refers to 
DM annihilation into neutrinos, while the right panels stands for annihilation into tau leptons. The 
elastic scattering cross section on protons (relevant for capture in the Sun) is fixed at the benchmark 
value of 10~ 41 cm 2 . The solid (red), dotted (blue) and dashed (black) lines refer to significance of 
2, 3 and 5 a. The two horizontal (pink) lines denote the exposures that can be reached by a 0.5 
Mton detector, like HK, in 1 (solid) and 10 (dotted) years. In this plot we assumed full efficiency 
for the detection of tau hadronic events (e T = 100%) and no misidentification of electron and muon 
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Figure 7. The same as in Fig. 6, for detection efficiency of tau hadronic events e T = 40% and 



misidentification of electron and muon events e™ 8 



4%. 



as an illustrative example, in Fig. 9 for a DM mass m x = 100 GeV and the benchmark value 
of a p . We notice that an efficiency of tau-events reconstruction of 50-60% would allow 
to a clear signal detection even without the need to reduce the level of misidentification. 
At the same time, an alternatively, a reduction of e™ ls and e™ ls at the 1% level would be 
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Figure 8. The same as in Fig. 6, for detection efficiency of tau hadronic events e T = 70% and 
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Figure 9. Dependence of the statistical significance <r on the misidentification parameters e e = e M 
and e T , for a DM mass m x = 100 GeV with cross section ct p = 10 _41 cm 2 and an exposure M^Ny = 
1 Mtonxyear. The solid (red), dotted (blue) and dashed (black) lines refer to <j = 2,3,5. The left 
panel stands for DM annihilation into neutrinos, while the right panels for annihilation into tau 
leptons. 



enough (for this benchmark case) to obtain a 5 a detection quite easily. 

In order to test the sensitivity of the downward-going tau signal not only to the DM 
mass, but also to its scattering cross section a p , we show in Fig. 10 the contours for 
q = 1.64 (which corresponds to a C.L. of 90%) in the plane a p vs. m x , for a detector with 
exposure M& et N y =1 Mtonxyear. The dotted lines represent the limits without considering 
misidentification, while the solid lines refer to e T =40% and e™ ls =e™ ls =4%. The left panel 
refers to spin independent interactions, while the right panel shows the case for a spin 
dependent cross section. These curves can alternatively be considered as the 90% C.L. 
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upper bounds deriving from the downward-going hadronic events on a p vs. m x . In the 
left panel of Fig. 10, beside our limits, we show the allowed regions obtained considering 
the DAMA [47-49], CoGeNT [50] and CRESST [51] positive results. The DAMA and 
CoGeNT regions are taken from Ref. [52] to which we refer for more detailed information. 
The DAMA region represents the domain where the likelihood function values differ more 
than 7.5 a from the null hypothesis (absence of modulation), while the CoGeNT region 
refers to the area where the likelihood function values differ more than 1.64 a. We also 
plot the constraints from XENON [53] and CDMS [54] experiments, as calculated in the 
statistical analysis of Ref. [55]. These limits are taken at 5 a and the XENON threshold 
has been set to 8 photoelectrons. In the right panel of Fig. 10, we report the DAMA 
region in the case of spin-dependent interaction [56]. We also report the limit on spin- 
dependent interactions from the SIMPLE experiment [57], see also Ref. [58] and Ref. [59] 
for discussions on critical points. Further weaker limits, not reported in Fig. 10, come 
from the PICASSO [60], COUPP [61] and KIMS [62] experiments. The recent analyses of 
direct detection annual modulation effects observed by DAMA [47-49] and CoGeNT [50] 
(and the excess of events reported by CRESST [51]) point toward a light DM candidate 
with a mass around 10 GeV and spin-independent scattering cross sections of the order of 
10 -42 cm 2 - 10 -40 cm 2 . For this type of particle, we would expect, for the case of direct 
annihilation into neutrinos, between 9 and 900 hadronic events and a detection with a 
statistical significance close to 5 a with a 10 years exposure on HK (5 Mtonxyr) in the 
case of the improved performance of Fig. 8. 

6 Conclusions 

In the context of indirect DM searches with neutrinos, the most common channel of in- 
vestigation is represented by upward-going muons or by contained //-like or e-like events 
produced by the charged-current conversion of the muon- or electron-neutrino fluxes pro- 
duced by DM annihilation in the Sun or in the Earth. These channels are very solid and 
do not suffer from large detection difficulties. In the most typical case of upward-going 
muons, the conversion area for the Vn — )■ fx process is represented by a large portion of the 
rock below the detector and experimental apparata posses very large detection efficiencies. 
For instance, SK has almost a 100% efficiency to detect through-going muons. The most 
important limit for this type of DM searches, instead, is represented by the large and 
v e atmospheric background that cast a shadow on the DM signal. 

In this paper we propose a new channel for DM searches at neutrino telescopes: 
downward-going hadronic tau events originated by the v T signal produced by DM an- 
nihilation in the Sun. This specific signal potentially represents a very good opportunity 
for DM, since the background of atmospheric downward-going v T is extremely reduced with 
respect to the upward-going z/„ case commonly considered. The intrinsic amount of v T in 
atmospheric neutrinos is very small as compared to v e and v„ components, while the v T 
component in the signal from DM annihilation (or decay) in the Sun is typically expected 
to be of the same order of their v e and counterparts. Moreover, in the downward-going 
direction, atmospheric do not have enough baseline to oscillated into v T . Therefore, a 
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m x [GeV] m x [GeV] 

Figure 10. Dependence of the statistical significance q on the DM scattering cross section on 
protons a p as a function of the DM mass m x , that can be derived by using the downward-going 
hadronic events with an exposure Md ct N y =1 Mtonxyear. The curves refer to ? = 1.64 and 
can alternatively be considered as the 90% C.L. upper bounds deriving from the downward-going 
hadronic events on a p vs. m x . The dotted lines represent the limits without considering misiden- 
tification, while the solid lines refer to e T =40% and e' e nis =e™ ls =4%. Left panel: our results for the 
spin independent case together with the allowed regions from DAMA (orange solid line for the case 
without channeling, orange dashed line for the channeling case), CoGeNT (dot-dashed red curve) 
and CRESST (cyan regions), and the limits from XENON 100 (green dashed line) and CDMS (gray 
dashed line) experiments. Right panel: our results for the spin-dependent case, together with the 
allowed regions from DAMA (upper for no-channeling, lower for channeling) and the limit from 
SIMPLE (gray dashed line). 

flux of tau neutrinos coming from the Sun when the star is above the horizon represents a 
signal with a very reduced background. The signal-to-background ratio in terms of fluxes 
is therefore much more favorable for down-going v^s then for up-going za/s and f e 's (or 
even v^s themselves). 

Additional sources of background are represented by v T produced in the solar corona 
of the Sun (which represent an irreducible background for our DM signal, since it comes 
from the same source, the Sun) and in v T arriving at the Earth from the galactic plane, 
mainly produced by oscillation of produced in cosmic-ray interactions with the inter- 
stellar medium. This source of background is, in principle, reducible by angular cuts on 
the galactic plane position. For detectors located in the northern hemisphere the angles 
cosQz > 0.5 are particularly favourable, because the galactic plane duty-factor is small 
for those angles. We have shown that all these sources of backgrounds are, by themselves, 
not posing relevant limits to the DM signal, which therefore represent a potentially viable 
novel possibility. 

The large advantage represented by the reduced background for the downward-going v T 
fluxes is, unfortunately, partly diminished by limits inherent in the difficulty to detect and 
properly identify tau neutrinos. Indeed, detectors specifically designed for the identification 
of tau events, like Opera at Gran Sasso Laboratory, are currently too small to collect a 
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statistically significant number of events. For this reason we have focused our analysis on 
Cherenkov detectors. In this case the hadronic tau events cannot be easily distinguished 
by the neutral current events, mostly coming from atmospheric v e and u^, and by the 
charged current e-like events. The possibility to correctly detect hadronic tau events is 
currently based on statistical methods and the percentage of misidentified events for a 
SK-like detector is of the order of several percent [36, 37]. As the atmospheric u e and 
are more abundant than the signal v T events, the misidentification influences the discovery 
potential of DM through the downward-going tau channel. 

Since for existing Cherenkov detectors, like SK, the number of hadronic tau events 
expected is small, we have focused our analysis on future Mton-size Cherenkov detectors, 
like Hyper-Kamiokande [11], UNO [63] or MEMPHYS [64], respectively in Japan, US and 
France. All these future detectors will be situated in the norther hemisphere and this could 
represent a great advantage to reduce the background of galactic neutrinos in the study of 
events from zenith angles cos#z > 0.5. We have not discussed in our analysis the IceCube 
detector, since v T reconstruction, at the low energies required to study the DM in the mass 
range of interest in our analysis (10 GeV - 1 TeV), may not be possible, but the possibility 
to implement statistical analyses is under study [65]. For the detection of high energy u T 
at IceCube see for instance Ref. [66] . 

For the Mton-size water Cherenkov detectors we have shown that the downward- 
going tau neutrinos signal has potentially good prospects, the main limitation being the 
level of misidentification of non-tau events, which need to be kept at the level of percent. 
For definiteness, we have studied two benchmark cases: DM directly annihilating into 
neutrinos, with equal amount of the three active flavors; dark matter annihilating into 
tau pairs. We showed that several tens of events per year (depending on the DM mass 
and annihilation/decay channel) are potentially collectible in a Mton-scale detector. Once 
the misidentification of non-tau events is taken under consideration, a 5 a significance 
discovery is potentially reachable for DM masses in the range from 20 to 300 GeV with a 
few years of exposure, and for a benchmark value of DM scattering cross section on protons 
(relevant for DM capture in the Sun) of a p = 10~ 41 cm" 2 . For light DM candidates with a 
mass around 10 GeV and spin-independent scattering cross sections of the order of 10~ 42 
cm 2 - 10~ 40 cm 2 , which are of special current interest due to the recent results in direct 
detection studies, we find that, for the case of direct annihilation into neutrinos, a detector 
like HK within 10 years (exposure of 5 Mtonxyr) could collect between 9 and 900 hadronic 
tau DM events and be sensitive to detection with a statistical significance close to 5 a, in 
the case of a 70% detection efficiency in the reconstruction of the tau hadronic events and 
of a 1% level of misidentification of non-tau events. 
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